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ABSTRACT 

A FORTRAN IV computer program was employed to conduct a 
statistical analysis of data collected during fleet anti- 
submarine warfare exercises. The object of the investigation 
was the identification of those variables which had greatest 
influence on a destroyer's ability to detect a submarine 
under certain conditions. 

The variables were treated as a random vector arising 
from one of two multivariate normal populations with common 
covariance matrix. An artificial regression relation was 
formulated to facilitate development of a linear discriminant 
function in a subset of those variables found to be of dominant 
importance. This latter subset was identified by examination 
of multiple correlation coefficients. 

The discriminant function was found to be seventy five 
per cent effective in classifying the experimental data 
correctly . 
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I. 



INTRODUCTION 



The sonar-equipped • destroyer is the primary surface 
vessel employed by the U.S. Navy in antisubmarine warfare. 
Detection of the submarine necessarily precedes any offen- 
sive action the destroyer can take against it. It was the 
objective of the author to identify those variables which 
had the greatest influence on the effectiveness of a 
destroyer in detecting a submerged submarine during fleet 
exercises. Consideration was limited to detection capabil- 
ities while using sonar equipment in the active mode. 

Reference 2 lists many of the large number of variables 
which influence the effectiveness of a destroyer in submarine 
detection. Certain of these variables can be regarded as 
control variables in the sense that the value of those vari- 
ables can be chosen and fixed by the destroyer commander with 
a certain degree of freedom. For example, the destroyer's 
speed may range from zero to thirty knots. Furthermore, a 
particular destroyer can be equipped with one of several types 
of sonar system, and can be manned with differing numbers of 
sonar technicians of various skill levels. These latter 
variables serve as examples of control variables whose value 
can be adjusted by agencies exogenous to the ship. Similar 
comments can be made concerning the submarine involved in an 
encounter with a destroyer. No attempt was made by the author 
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to consider any variable as a control variable; all were 
regarded as completely random variables. 

Based on personal experience in the Destroyer Force, the 
author chose a set of twenty-two variables whose examination 
was expected to be most fruitful for the purpose stated. 

Four groups of variables were considered; environmental 
factors, destroyer characteristics, submarine characteristics, 
and tactical factors. Table I lists the variables included 
in each group. Qualitative variables, such as destroyer hull 
type were assigned numerical codes. All values were recorded 
at the time detection was achieved, or - in those cases where 
detection did not occur - when the submarine had approached 
to a specified range. This range was specified by the U.S. 
Navy Antisubmarine Warfare Data Center (NADAC) as a standard 
for each sonar system. The standards are listed in Ref. 2. 



Table I . Variable?; by Group 



Environmental 


Destroyer 


Submarine 


Tactical 


Geographic 

location 


Hull type 


Hull type 


Relative speed 

between 

opponents 


Present weather 


Hull number 


Hull number 


Distance between 
opponents 


Sonic layer 
depth 


Sonar type 


Depth 


Aspect angle of 
submarine 


Wind speed 


Speed 


Speed 


Relative bearing 
from destroyer 
to submarine 


Sea state 


Self noise 






Surface turbu- 
lence 


Technician com- 
plement 






Bathythermo- 
graph profile 
type 








Water depth 
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The data examined was supplied by NADAC on request, having 
been collected during sixteen fleet exercises conducted in 
the Atlantic and Pacific Oceans during the period 1962 to 
1966, inclusive. The destroyer types considered included DLG, 
DL, DDG, DDR, DD, DEG and DE. Each destroyer was equipped 
with one of the following AN/SQS sonar systems: 23, 26, 29, 

or 30. Submarine types encountered included all types in the 
U.S. fleet during this period, ranging from the World War II 
vintage fleet submarine to the nuclear-powered ballistic 
missile submarine. 

Each observation of the twenty-two variables was regarded 
as a random, vector, x, arising from one of two multivariate 
normal populations with common covaariance matrix: N(mj^, V), 
or N (m 2 , V) . The former distribution representing an encounter 
which resulted in detection of the submarine by the destroyer, 
and the latter distribution applying otherwise. 

Maximum likelihood estimators were used for the mean 
vectors, covariance matrix, and all other parameters employed 
in the analysis. 

a FORTRAN IV computer program was prepared for analysis 
of the data and was utilized with the U.S. Naval Postgraduate 
School IBM-360 computer. A discussion of the program, a flow 
chart, and the program appear below. 

In all, observations from ninety two encounters were 
included in the analysis. Twenty-seven encounters resulted in 
submarine detection by the destroyer, and sixty-five did not. 
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II. DISCUSSION 



• A. PRINCIPAL STEPS OF THE ANALYSIS 

The following discussion is concarned with the principal 
concepts employed in the analysis , without regard to the 
computer programming techniques employed. A full discussion 
of the program appears later. 

The linear discriminant function 

x'j v“^ (m^ - m 2 ) - y(m^ + m2)' v“^ (m^ " ® 

is related to the artificial regression proposed by Fisher 
in Ref. 5 

Vj = k' “ £) • (2) 

In the foregoing, is the observation vector with j = 1 , 

2 , ... , n^ + n2 ; there are n^^ observations in the first 
population with parameters V) and n2 observations in 

the second population. The artificial dependent variable, 
y., takes the value n_/(n^ + n_) if x. is from population 
one, and the value -n^^/Cn^ + n2) otherwise. The vector c 
is a constant vector with each element being the weighted 
grand mean of observations, i.e., (n^^ x^^ + n2 X2)/(nj^ + ^ 1 ^ 2 ^ ' 

Equation ( 2 ) yields artificial regression coefficients, 
b, which are proportional to the discriminant function coef- 
ficients, V ^ (ra^ - m2) / of equation ( 1 ) . (This may be 
verified by solution of the normal equations in the components 
of b.) 
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Based on the above relationship between the discriminant 
and regression functions, equation (2) was utilized to 
determine a subset of variables which could be regarded as 
being of principal importance in the discriminant function. 

The technique employed is discussed in Ref. 7. A multiple 
correlation coefficient, Rj.!' calculated from equation 

(1) using the full set of twenty-two independent variables. 
Then, each independent variable, x. , was eliminated in 
iterative fashion and a corresponding multiple correlation 
coefficient, 1 ^ 21 ' calculated based on the remaining twenty 

one independent variables. The significance of each of the 
twenty-two variables in the regression was then tested 
according to the following scheme. If the following relation- 
ship was satisfied, then the deleted variable, , was discard 
ed as contributing little to the regression: 

q(Rll - " ^11^ ^Pl " ^ 

In this relationship p^ was the number of variables in the 
original regression (22) , P 2 was the number in the reduced 
set (21) , F was the F variate with (p^^ - P 2 ) and q degrees 
of freedom, respectively; q had the value (n^ + n 2 ) ~ 

+ 1 ) . 

The above technique revealed that all but six of the 
original variables contributed little to the regression, and - 
therefore - to the discriminant function. The value of R^^ 
was found to be 0.646 for the original variable set and had 
the value 0.528 for the six significant variables. The latter 
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value was calculated from a regression based on these six 
variables alone. That is, with a seventy-three per cent 
reduction in the number of independent variables, the corre- 
lation coefficient was reduced by only eighteen per cent. 

Table II lists the six significant variables in order of 
importance, the value of the correlation coefficient for the 
set of twenty-one variables remaining upon deletion of each 
of the significant variables , and the reduction in the multiple 
correlation coefficient which resulted upon deletion of the 
variable . 

With the six significant variables thus identified, a 
linear discriminant function was constructed, and its properties 
tested. The specific form of equation (1) used at this point 
was 



x’ d - ik = 0, 

-D - 2 



(3) 



,th 



where x. is the i 6X1 observation vector, d is the vector 

-j J 

of discriminant coefficients and k is a constant. The relation 
between equations (1) and (3) is given by 



d = V ^ (m^ - 1 ^ 2 ) 

k = (m^ 512^ ' ^ ^ — 2 ^ 



With the value of all parameters calculated, equation (3) 
took the follov7ing form: 



-O.ieSXj^ + 0.015x2 + 0 . 411 x 2 - 0.020x^ + O.SOGx^ 

- O.OOSx^ > 8.28 . 
b — 



(4) 
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Table II. Significant Variables. 



i 


X . 
1 


R„. 

2r 


Per cent Reduction 
in Correlation 


5 


Type of weather 


0.565 


12.5 


2 


Submarine hull number 


0.604 


6.5 


6 


Sonic layer depth 


0.606 


6.2 


1 


Geographic location 


0.608 


5.9 


3 


Submarine hull type 


0.612 


5.3 


4 


Submarine depth 


0.617 


4.5 



An observation would be classified as arising from population 
one if the relation was satisfied, and as arising from popula- 
tion two otherwise. With the formulation of equation (4), the 
observations were then to be classified using the equation, 
and the results achieved to be considered a measure of the 
validity of the discriminant function. 

Anderson discusses in Ref. 1 the probability of mis- 
classif ication associated with a linear discriminant function. 
Under the assumption that the cost of an error of misclass- 
ification is equal for each of the two possible errors, the 
probability of misclassification, P (M) , is given by 

P (M) = P (Y > /a/2' ) 

where Y is a standard normal variate and a takes the value of 
(m^^ - 1 ^ 2 ) V ^ (m^ “ 'the above equation, it was 

found that the value of P (M) associated with equation (4) is 
0.18. 

Upon classifying the ninety-two observations of the six 
significant variables - using equation (4) - it was found that 
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eight (30%) of the twenty-seven cases where detection occurred 
were misclassified as cases of non-detection. Similarly, 
fourteen (22%) of the sixty-five cases of non-detection were 
misclassified. Thus, the overall error of misclassif ication 
was twenty-five per cent. 

In view of the small sample sizes available, particularly 
in the case of population one (twenty-seven observations) , the 
difference between the theoretical error of misclassif ication 
and the observed value was considered acceptable. The fact 
that the error of misclassif ication in each population 
decreased with sample size was viewed as an encouraging sign. 

Consequently, it was concluded that the six variables 
cited above were, in fact, of dominant importance in destroyer - 
submarine encounters. 

It was noted with interest that the character of these 
variables is such that the destroyer commander has virtually 
no m.eans available to influence them to his advantage. There 
seems to be little direct action that could have been taken 
by the destroyer commander to increase significantly the like- 
lihood of detecting his submerged opponent. 

On the other hand, the submarine commander enjoys a distinct 
advantage. If his depth and location are chosen judiciously, 
and his submarine is of the preferred design, he may proceed 
with confidence of not being detected by his surface adversary. 

B. THE COMPUTER PROGRAM 

The computer program employed in the analysis contained 
four principal steps. First, the data were read into storage. 
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constants were defined, and the data were transformed into a 
form to facilitate subsequent calculations. Second, the 
multiple correlation coefficient for the twenty-two independent 
variables was calculated. Third, using the scheme described 
above, the significant variables were identified. Finally, the 
discriminant function was formulated and tested. The discus- 
sion below is concerned with the mechanics and logic of the 
entire program. A flow chart of the program is included as 
Appendix A, and is followed by a listing of the program itself. 

The following subprograms which were prepared by the 
computer manufacturer formed a substantial part of the program: 
CORRE, ORDER, MINV, MULTR, DMATX . A full discussion of these 
subprograms may be found in Ref. 6. They are discussed below 
briefly where necessary to maintain continuity. 

After storage requirements and constants are defined in the 
main program, control is passed to the subroutine DATA. Here, 
a DO loop is used to read the raw data into storage as a 
three-dimensional array X(I,J,K). (The indices of X are 
observation number, variable number, and population number 
respectively.) As each card of data is- read, the values of 
four variables are calculated and the data are transformed to 
a one dimensional array XD. The latter step is required by 
the svibsequent use of the manufacturer-supplied subprograms 
mention above . 

The first value calculated is that of the dependent var- 
iable y^ which is stored in X(I,23,K). Next the relative 
bearing from destroyer to submarine is converted to 
radian measure using the simple subprogram REL. With destroyer 
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course and submarine true bearing as input variables, REL 
calculates the relative bearing as the output variable RB. 

This value is then converted to radian measure and stored in 
X(I,12»K) . Then, relative speed is computed from destroyer 
speed, submarine speed, and the angle between the two vessels' 
courses. This value is stored in X(I,5,K). The value of 
sonic layer depth is stored in X(I,22,K) based on the fleet 
exercise from which the observation was drawn. Finally, the 
destroyer's manning level in the sonar technician (ST) rating 
is used to calculate a parameter reflecting the rate distri- 
bution. This is accomplished by calculating a sum of v;eighted 
values and storing the result in X(I,13,K). For each destroyer, 
the number of chief petty officers is multiplied by six, first 
class petty officers by five, and so on, until a product is 
formed for all skill levels in the ST rating. The sum of these 
products is then taken as the measure of manning level for that 
destroyer . 

With the data transformation completed, control returns 
to the main program and the multiple linear regression 
equations are formulated. The first subprogram called is 
CORRE. Given as input variables the number of observations 
NS, the number of variables M, and the one dimensional array 
of observations, XD - CORRE computes the vector of means 
CMEAN, the vector of standard deviations PS, the MxM matrix of 
sums of cross-products of deviations from means D, and the 
upper triangular portion of the MxM symmetric matrix of 
correlation coefficients R. 
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After the subscripts of the variables are stored in the 
array ISAVE in the order in which the variables are arranged 
in the array XD, control is passed to siibroutine ORDER. 

Given the nuinber of variables, the matrix of correlation 
coefficients, the subscript number of the dependent variable, 
and the array ISAVE, subroutine ORDER computes the 22x22 matrix 
of correlation coefficients among independent variables 
and the vector of correlation coefficients between the dependent 
variable and the independent variables R 22 * 

Next, the subroutine MINV is called to calculate the 
inverse of This matrix plus the vector of means, the 

vector of standard deviations and the diagonal elements of 
the matrix D are then passed to the subroutine MULTR. This 
subprogram computes the vector of regression coefficients W 
and the multiple correlation coefficient ANS(2). (The vector 
W and ANS(2) correspond to the vector b of equation (2) and 
the variable R^^^ of the relation (2ci) discussed earlier.) 

A DO loop is then used to eliminate each independent 
variable from the regression, and twenty-two multiple corela- 
tion coefficients for each combination of twenty-one variables 
is computed. The test statistic of relation (2a) is computed 
for each variable x^ and stored in C(J) . As each new 
multiple correlation coefficient is calculated, the value of 
the test statistic is stored temporarily in TEM3 for each 
test. This value is then compared to the value of the F(l,69) 
variate at the 5% significance level. If the variable with 
subscript J is found to be of little significance, then a zero 
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is stored in array MK, otherwise MK(J) will contain the value 
of the subscript. 

On completion of the above, the values of the test sta- 
tistic are arranged in order of magnitude in the array C, and 
deleted from the array ISAVE. 

At this point the number of significant variables is 
checked. This step is necessitated by the fact that the array 
TST contains values of the F(f^,f 2 ) statistic with f^^ = 1 and 
6 _< f 2 ^ 14. The statistic is used later in the program and 
the results of the program would be invalid if the niamber of 
significant variables were not within the range given for f 2 . 
(The range chosen was based on preliminary trials with the 
program. It was expected that there would be six significant 
variables . ) 

The subprograms ORDER, MINV and MULTR are employed as before 
with the regression equations being based on the significant 
variables only. Deletion of the indices of the insignificant 
variables from the array ISAVE results in the appropriate 
adjustments in the subprograms. The test statistic TEM3 is 
compared to the appropriate F variate TST (IT), and the index 
is incremented. If the test is not satisfied, i.e., if the 
discarded variables have significance as a group (an unexpected 
result) , a value of one is stored in ICK and this will result 
in a change of the flow of the program later. 

A DO loop follows the above and rearranges the order of 
the variables in the array XD in preparation for calculation 
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of the discriminant function. Control then passes to the 
sxibroutine DMATX. 

Given the nixmber of (significant) variables L, the number 
of populations LX, the number of observations in each popula- 
tion N, and the array of data XD - DMATX calculates the vector 
of means XBAR and the common covariance matrix D. The inverse 
of D is then calculated by subroutine MINV and control passes 
to subroutine DISC. This latter subroutine is a simple 
FORTRAN IV equivalent of equation (1) . 

Next, each observation is classified using the discriminant 
function just formulated. The value of the function is tested 
against zero and the observation classified accordingly. For 
each population, the value of KOUNT is incremented each time 
an observation is improperly classified. 

After all observations have been classified, the value of 
ICK is checked. As discussed above, if ICK has value one, the 
previously discarded variables have group significance. If 
this is the case, a warning to that effect is printed, the 
index of the most significant deleted variable is added to the 
array MK, control returns to statement number 45, the new 
group of significant variables are tested and a new discriminant 
function is found. If ICK has a value of zero, the discrim- 
inant function coefficients are printed and the program 
terminates . 
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III. CONCLUSION 



Analysis of data collected during fleet exercises conducted 
during the period 1962 to 1966 inclusive strongly suggests 
that the variables listed in Table II had a dominant influence 
on destroyer - submarine encounters. Specifically, the 
successful detection of a submarine by a destroyer utilizing 
a sonar of the type specified in the active mode v/as influenced 
more strongly by these six variables than by the remaining 
sixteen considered. 

The discriminant function, equation (4) , was 75% effective 
in classifying encounters and should be useful as a predictive 
tool in assessing the likelihood of detection under knov;n 
environmental conditions and against, known submarine types. 

Consideration of the variables found to be significant 
indicates that the tactical advantage was with the submarine 
commander during encounters of the type considered. 
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FLOW CHART OF COMPUTER PROGRAM 
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COMPUTER PROGRAM 



DIMeNSION X{ 65 ,23 »2 ) , X3AR(44 ) , 0(529), CMEAM(23), 
/XD(2116) ,PS(23) ,M(2) ,R(529) ,3(23) ,W(25) ,T(23) , 

/ ISAVF( 2 3 ) , PI 1 ( ^^84 ) , R 1 2 ( 22 ) , S ( 2 3 ) , AN S ( 1 C ) , MK ( 22 ) , 
/IDX(22) ,TST(10) ,C(23) 

RCAD(5,1) LX,M,N 

ICK=C 

IT=0 

KI=P-1 

NS=N(1 ) +N( 2 ) 

DATA T ST /2. 36, 2. 39, 2. 42,2.45,2.49,2.55,2.60,2.68, 
/2.8C,2.Q5/ 

C 

C TRA^■SFGRM PAV! DATA FOR MULTIPLE REGRESSION ANALYSIS 
C 

CALL DATA( M,N,MS,X,XD) 

CALL CORRE (NS, M, 1 , XD,CMEAN, PS ,0 ,R ,B ,W,T ) 

DO 1C 1=1, K I 
10 I SAVE( I ) =I 

CALL CRt'LR (N,R,M,KI , I SA VF , R 1 1 , R 1 2 ) 

CALL MINV( PI 1 , KI , DT ,W, T ) 

CALL MULTP (NS,KI , C ME AN , PS , B , R li , R 1 2 , 1 S A VE , W , S W , T , ANS ) 
WRITE (6 ,8C C ) 

WRITF(C,5C 1 ) ( W( I A) , I A=1 , KI ) 

WRITE! 6, 81 C ) 

WRITE(6,5C1) ANS(2) 

TFMl = (AMSt 2) )’^’^2. 

DN=( l.-Ttf 1 ) /69. 

KD=KI-1 
DO 30 J=1,KI 
L = 0 

DO 2C JI=1 ,KI 
IF(JI.EO.J) GO TO 20 
L=L + 1 

ISAVE (L)=JI 
20 CONTINUE 

CALL ORDER (M,R, M,KD, I S A VE , R 1 1 , R 1 2 ) 

CALL MI NV( R1 I , KD ,DT, W, r ) 

CALL MULTR (NS, KD,CMEAN ,PS,B,R11 ,R1 2 , I S AVE , W, SW , T , AMS ) 
WRITE. (6,820) J 
WRITE(6,5C1) ANS(2) 

TEM.2=(ANS( 2) )«-*2. 

TLM3=(TEM1-TEM2 )/DN 
C( J) =TEM3 

IF(TEM3.GT.3.98) GO TO 27 
25 MK(J)=0 
C 

C MK(J)=0 IMPLIES THAT X(J) DOES MOT CONTRIBUTE 
C SIGNIFICANTLY TO THE CORRELATION COEFFICIENT. 

C 

GO TO 3C 
27 MK(J)=J 
30 CONTINUE 

DO 4C KJ=1 , 22 
STR=C( KJ) 

IND=KJ 

DC 35 KL=1,22 
IF(C(KL).LE.STR) GO TO 35 
STP=C( KL) 

IND=KL 
35 CONTINUE 

I D X ( K J ) = I N D 
40 C( IND)=C. 

C 

C IDX CONTAINS THE INDICES OF THE ' ' I NDEPENDENT • • 
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C VARIABLES ORDERED WITH RESPECT TO SIGNIFICANCE OF 
C THE VARIABLE’S CONTRIBUTION TO THE CORRELATION 
C COEFFICIENT 
C 

A5 L = 0 

DO 60 JJ = 1 ,KI 
IF(MK( JJ ), EQ.O ) GO TO 60 
L = L+1 

ISAVE( L)=MK( JJ) 

60 CrhTINUE 

IF( ( L, LT. 6) .OR. IL.GT, 15) ) GO TO 1010 

WRITE! 6, 83 0) 

WRITE (6, 2) ( ISAVE( IV) , I V=li L ) 

CALL ORDEP (M,R,M,L, I S A VE , R 1 1 , R 1 2 ) 

CALL MINV(ril,L,DT,W,T) 

CALL N’JLTRINS, L , CMEAM, PSi B, Rll , R12, ISAVE , W, SW,T, ANS ) 
WRITE(t,8AC) 

WRITE(6,501) ANS(2) 

TEN2=(ANS( 2) >’'='^2. 

TFM3=( TER1-TEM2 ) / ( DN*- ( K I -L ) ) 

IT=IT+1 

IF(TEM3.GT.TST ( IT) ) ICK=1 

DISCRININAMT FUNCTION IS FORMED BASED ON THE 
SIGNIFICANT • • I ND EPENOENT ’ ’ VARIABLES. 

LK = 0 

DO eo K=-l,2 
NK=N(K) 

DO 70 JK=1,NK 
LK=--LK+1 
N2=LK-NK 
DO 70 IE=1,L 
JR=ISAVE(T E) 

N2 = E 2+NK 

70 XD(N2)=X( JK, JR,K) 

80 LK=N2 

CALL DM ATX ( LX i L t N , XDf XB AR t D t T ) 

WRITE (6, 99 5) 

DO 77 I 1=1 , L 
JX=I I+L 

77 WRITE! 6 ,99*=) X B AK ! I I ) , XBAR ! J X ) 

CALL MI NV! D, L, DT, T, C) 

CALL DI SC! XBAR ,L,D, T,C ) 

INPUT DATA IS CLASSIFIED ACCORDING TO THE CURRENT 
DISCRIMINANT FUNCTION 

DO 81 KN=1 ,2 
KOUNT=0 
N1=N!KN) 

DO 85 I L=1 , N1 
SUM=0. 

IND=C 

DC 79 JL=1 , L 
IG=I SAVE! JL) 

IND= IND+1 

7 9 SUM=SUM + C ! riD)«X! IL,IG,KN) 

IF!SUM-C!l.+1 ) ) A, 5, 5 
A NM=2 

GC TO 88 
5 NN=1 

88 I F !NN.L,E.KN) KOUNT = KOUNT + 1 
85 CCNTINUE 

81 WRITE !6,6C 1 ) KOUNT,KN 
IF ! ICK. EC. C ) GO TO ICCO 

90 WPITE!b,91 ) 

IND= IDX !L+1 ) 

MK ! I MD ) = I N D 
ICK = 0 
GC TO ^5 
1000 LK=L+1 
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WRITE( 6, 500 ) 

WRITE(6,501) ( C( JS) , JS=1,LK> 

STOP 

1010 WRITE( b, 10 15 » L 

1015 FORMAT ( 'WARMING* , ' L=* t I 10 ) 

STOP 

1 FORMAT 1514) 

2 FORMAT { lOX ,2014) 

91 FORMAT ( lOX , *01 SCARDED VARIABLES HAVE GROUP 
/•SIGNIFICANCE' ) 

500 FDRMAT(25X , 'DISCRIMINANT FUNCTION COEFFICIENTS' ) 

501 FORMAT ( 8E1 5. 3) 

601 FORMAT ( IHC , I 35 , 4X , • M I SCL ASS I F I CAT I ONS IN GROUP', 16) 
800 FORMAT ( IHC , lOX ,' STANOA RD REGRESSION COEFFICIENTS') 

810 FORMAT ( IhC , 'MULTIPLE CORRELATION COEFFICIENT') 

820 FORMAT ( IHC ,' MULTI PLE CORRELATION COEFFICIENF FOR ', 
/•GROUP • , 14 ) 

830 FORMAT! IHO, • INDICES OF VARIABLES MAKING SIGNIFICANT ', 
/• CCrnR I BUT ION TO THE CORRELATION COEFFICIENT') 

340 FORMAT! IHC, 'CORRELATION COEFFICIENT FOR SIGNIFICANT ', 
/ • VAR lABLES ' ) 

995 FORMAT !28X ,' MEAN VECTOR ONE ' , 25 X , • ME AN VECTOR TWO') 

999 FORMAT ! 2E4 0,3) 

END 



SUBROUTINE REL! PHI 1,PHI2,PHI ) 
PH 1= PH I 1-PHI2 

IF!PHI,LT.O.O) PHI=PHI+360. 
RETURN 
END 



SUBPOUTINF DISC ! XBAR , M , D , P S , C ) 
DIMENSION XBAR ! 1 ) , PS! 1 ) ,C ! 1 ) ,0! 1 ) 
SUM=C. 

DO 76 J=l, M 
MC=J-*M 

76 PSIJ) =XBAR ! J) -XBAR!MC ) 

IJ = C 

DO 77 1=1, M 
C! I )=0. 

DO 77 L=l, M 
IJ=IJ+1 

77 C!I )=C! I)+PS!L)^D! IJ) 

DO 76 K=l, M 

MC=K+M 

78 SUM=SUM+ !X BAR ! K ) 4-XBAR! MC) )^C! K ) 

C!M+1.)= !SUM/2.) 

RETURN 

END 



SUBRCUT INF M I N V ! A , N , D , L , M ) 

DIMENSION A! 1 ) ,L! 1 ) ,M! 1 ) 

D=1.0 

NK=-N 

DO 80 K=1 , N 
NK=NK+N 
L !K)=K 
M! K) =K 
KK=NK+K 
B IGA=A! KK ) 

DO PC J=K,M 
IZ=N=«=! J-1 ) 

DO 20 I=K,N 
IJ=I Z+ I 

10 IF! ABS!BIGA)- ABS!A!IJ))) 15,20,20 

15 BIGA=A! IJ ) 

L!K)=I 
M! K) =J 

20 CONTINUE 
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J=L( K) 

IF(J-K) 35,35,25 
25 KI=K-N 

DO 30 1=1, N 
K I=KI + N 
HGLD=-A(KI ) 

JI=KI-K+ J 
A( KI ) = A( JI ) 

30 A( JI) =HCLO 
35 I=M(K) 

1F(I-KI 45,45,38 
38 JP = N=«^( I -1 ) 

DO 4C J=1 , N 
JK=NK+ J 
J I=JP+J 
HOLD=-A( JK) 

A( JK)=A( JI ) 

40 A(JI) =FULr 

45 IF (BIGA I 48,46,48 

46 n=c.e 
WRITE(6,4C0» 

400 FORMAT (4CX SI NGULAR INPUT TO MINVM 
RETURN 

48 DO 55 1=1, N 

IF(I-K) 50,55,50 
50 1K=KK+I 

A( IK) = A( IK )/(-BIGA) 

55 CONTINUE 

DO 6 5 1 = 1, N 
I K=NK+I 
HOLD=A( IK) 

IJ=I-M 
DO 65 J=1,N 
1J=IJ+N • 

IF(I-K) 6C,b5,60 
60 IF(J-K) 62,65,62 
62 KJ=1J-1+K 

A( IJ ) = )iOLDvA( KJ )+A( IJ ) 

6 5 CCiNlINUE 
KJ=K-N 
DO 75 J=1,N 
K J = KJ+ N 

IF(J-K) 7C,75,70 
70 A(KJ )=A(KJ )/BIGA 
75 CONTINUE 
D = D4-BIGA 
A(KK )=1,C/BIGA 
80 CLNllNUE 
K = N • 

100 K=(K-1 ) 

IF(K) 150,150,105 
105 1=L(K) 

IF(I-K) 120,120,108 
108 JO=N*{K-n 
JR = N>^( I -1 ) 

DO no J = 1,N 
JK=JC+ J 
HOLD=A( JK) 

JI=JR-» J 
A( JK)=-A( J I ) 
no A( JI ) =UCLD 
120 J=M(K) 

IF(J-K) 100,100,125 
125 KI=K-N 

DC 130 I=1,N 
KI=K I + N 
HOLD=A( KI ) 

J I=K I-K+J 
A(KI ) = -A( J I ) 

130 A( JI ) =OCLD 
GO TO ICO 
150 RETURN 
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END 



SUBRCUTINF DHATX ( K , H , N , X , XBAR , D, CMEAN ) 
DIHFNSILN N( 1 ) ,X( 1 I ,XBAR( 1 ) , 0( 1 ) ,CMEAN( 1 ) 
MN=MfrM 

00 100 1=1, MM 

100 D( I )=0«0 
NA=0 
L=0 
LM=0 

DO 160 NG=1,K 
N1=N(NG ) 

FN = N1 

DO 130 J = 1 ,M 
LM=LM+1 
XRAR(LM)=0.0 
DU 120 1=1, Ml 
L = L + 1 

120 XDAR(LM)=XRAR(LM)+X(L) 

130 XBAR (Ll-M =XPAR( /FN 
LHFAN=LM-M 
DO 130 1=1 ,N1 
LL=NA+ I -ra 
DO 140 J=1,M 
LL = l.L+Ml 
M2 = l MFAN+J 

140 CMEAN(J)=X(LL)-XBAR(N2) 

LL = 0 

DC 130 J=1,M 
DO 130 JJ=1,M 
LL=LL+1 

150 D(LL) = 0(LU <-CMEAN( J )>''CMEAN( JJ ) 

160 N4 = n^+N1^-M 
LI =-K 

DO 170 I=1,K 
170 LL=LL+N(I) 

FF-LL 

DC 180 1=1, MM 
180 D( I )=D( I )/FN 
RETURN 
END 



SUBROUT I NE D AT A ( M , N , NS , X , XD ) 

D 1 MENS I ON N ( 2 ) , X ( 6 3 , 2 3 , 2 ) , SOD ( 8 , 66 ) , S LD ( 1 6 ) , XD ( 2 1 1 6 ) , 
/C(2) 

READ (3,200) SOD 

SOB IS. AN II RUT MATRIX CONTAINING DATA ^ELATING 
TO MANNING STRENGTH IN THE ST RATING FOR EACH 
DESTROYER PARTICIPATING IN THE EXERCISES . 



THE INPUT MATRIX SLD CONTAINS SONIC LAYER DEPTH 
DATA FUR EACH EXERCISE. 

READ (3,100) SLD 

C( 1)=FLCAT(N(2) ) /FLOAT(NS) 

C(2)=-FL0AT(N( 1 ) )/ FLOAT (NS) 



CONVERT 180 DEGREES TO RADIAN MEASURE 

DEN= 180. /3. 141593 
L = 0 

DO 68 K=l,2 
NK=N(K ) 

DO 61 1=1, NK 
L = L+1 
N2=L-NS 

READ(5,102) (X ( I , J, K) , J=1 , M) 

X(1 ,23,K) = C(K) 
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CIE*i 






c 

C CONVERT GEOGRAPHIC DATA TO RELATIVE DATA. 

C 

CALLRtK X< I , 12,K) ,X( I ,5,K) ,RB ) 

X( I , 12, K) = P3/D£N 

CALL PfiLiX ( I ,13 ,K) ,X( I , 5 , K ) , OUT ) 

OUT=GUT/DE N 

XI I,5.K )= < X( I,6,K)>;'=f-2.+X( I, l4,K)^'-*2.-2.-~X(I ,6,K)* 
/X( I ,14,K)<^C0S( OUT) 

DO 63 KK=1 , lo 
Y=FLOAT ( KK ) 

IFIXd , 1,K).EQ.Y) GO TO 62 
GO TO 63 

62 X(I ,22 ,K) = SLD(KK) 

63 CCMTINUE 

DO 64 t<D=l,150 

lEMXn, l,K),EQ,SOD(l,KD)).AND.(X{I,3,K).EQ. SOD( 2,KD) 
/ ) ) GO TO 7 

64 CCNT INUE 
7 P=6. 

X( I , 13, K) = 0. 

DO 74 KX=3,3 

X( I ,]3,K) = xn ,13,K)+P=i'S0U(KX,K0) 

74 P=P-1. 

DO 61 J J=1 , M 
N2=N2+ MS 

61 XD(rJ2)^X(I,JJ,K) 

68 L=NK 

100 FORMAT ( 7F1 0. 0) 

102 FORMAT (3F1C.0, lOX, 3F 1 C . G , / , 7F 1 0. 0 , / , 7F 1 0. C , / , F 1 0. 0 , 
/lCX,2r IC.C ) 

200 FCRMAT(F4. 0,F6.0,6F10.0) 

RETURN 

END 



SUBROUTINE CORPE ( N , M , I 0 , X , XBAR , STD , R X , P , D , D , T ) 

D1 MENS I ON X ( 1 ) ,XBAR( 1 ) , ST D ( 1 ) , P X ( 1 ) , R ( 1 ) , b ( 1) , 0 ( 1 ) , 
/T< 1) 

DC 100 J=1,M 
B( J) =C. C 
100 T(J)=C. 0 

K= ( f-R" M + M )/2 
DO 102 I =i ,K 
102 R( I )^0.0 
FN = N 
L=0 

IF(in) 105, 127, 105 
105 DO ICB J=1,M 
DO 107 I=1,N 
L = L+1 

107 T( J) =T(J )+X( L) 

XBAR(J)=T(J) 

108 T ( J )=T ( J )/FN 

DO 115 I=1,N 

JK=0 

L=I-N 

DO lie J=1,M 
L = L + N 

D< J)=X(L)-T{J) 

110 B< J) =0 ( J) + D( J) 

DO 115 J=1,M 
DO 115 K=1 , J 
JK=JK4 1 

115 R(JK ) = k( JK )4-D( J)vD(K) 

GO TO 205 

127 IF(N-M) 13C, 130, 135 

130 KK=N ( 

GO TO 137 
135 KK=M 

137 DC 140 I=1,KK 
CALL DATA (M,D) 
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DO 14C J=1,M 
T(J)=T( J)+D( J) 

L = L + 1 

140 RX(L)=J(J) 

FKK=KK 

DO 1^0 J=1 ,M 
XBAk ( J ) =T( J) 

150 T(J)=T( J)/FKK 
L=0 

DO leo 1=1 ,KK 

JK = 0 

DO 17C J=1,M 
L = L+1 

170 D(J)=FX(L)-T(J ) 

DO leo J=1 ,M 
B( J ) = B( J > + D( J ) 

DO 18C K=1,J 
JK = JK+ 1 

180 R( JK ) = ■• ( JK )+D( J )=^^D(K) 

IF(O-KK) 2C5, 205, 185 

185 KK=0-KK 

DO 200 1=1, KK 

JK=0 

CALL DATA (M,D) 

DC 190 J=1,M 

XBAKI J) =XBAR ( J ) +D( J ) 

D( J )=D( J )-l ( J I 
190 B( J)=B( J)+D( J) 

DO 200 J = 1 ,M 
DO 2C0 K=1,J 
JK=JK+1 

200 R{ JK ) = H' ( JK ) +0( J )*0 ( K) 

205 JK=C 

DO 2 1C J=l , M 
XBAM J ) =XBAR ( J ) /FN 
DO 2 1C K = 1 , J 
JK=JK+1 

210 R(JK ) = R(JK)-D( J)’^B(K)/FN 
JK=0 

DO 220 J=1,M 
JK=JK+ J 

220 STD(J)= SQRT( ABS(R(JKI)) 

DO 230 J=1,M 
DO 230 K=J,M 
JK = J + ( K>: K-K ) /2 
L=N'!' ( J-1 )+ K 
kX(L ) = P ( JK ) 

L=fi^= (K-1 ) + J 
RX(L).= P ( JK) 

IF ( STD( J STD( K ) ) 225, 222, 225 
222 k(JK)=0.0 
GO TO 230 

225 k.(JK ) = k(JK)/(STD(J)*STD(K) ) • 

230 CCOTIMUE 

rN=sokT(rc-i,o) 

DC 2^0 J=1,M 
240 STD( J)=STD( J)/FN 
L = -M 

DO 250 1=1 ,M 
L=L + M+1 

2 50 B( 1 ) =kX( L) 

RETURN 

END 



SUBROUTINE ORDER ( M , R , NDE P , K , I S AVE , RX , RY ) 
DIMENSION R( 1 ) , ISAVE( 1 ) ,RX( 1 ) , RY( 1 ) 

MM = 0 

DO 130 J = 1 ,K 
L2 = I SAVE ( J ) 

IF(NDEP-L2) 122, 123, 123 
122 L = NDCP+ (L2'^L2-L2 )/2 
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GO TO 125 

123 L=L2+(NDEPXNDEP-NDEP)/2 
125 RY(J)=R(L) 

DO 13C 1=] ,K 
L1=ISAVE(I ) 

IF(L1-L2) 127, 123, 128 

127 L=Ll+( L2VL2-L2 ) /2 
GO TO 129 

128 L=L2+( L1*L1-L1 )/2 

129 MM=MM+1 

130 RX(MM»=R(L) 

ISAVE( K+l) =NDEP 
RETURN 

END 



SUBROUTINE MULTR ( N, K , XBAR , STD, D, RX , RY , I S AVE , B , SB , T , 
/ANS ) 

DIMENSION XB AR< 1 ) , STD( 1 ) ,D( 1 ) ,RX( 1 ) ,RY{ 1 ) , ISAVEI 1 ) , 
/T ( 1 ) , ANSI 1 ) , B( 1 ) , SB I 1 ) 

MM=K+1 

DU 100 J=1 ,K 
100 B(J)=0.( 

DC 1 10 J = 1 ,K 
L1=K*( J-1 ) 

DO lie I=1,K 
L=L1+I 

110 B( J) =B( J) + KY( I I’i'RXI L) 

RM=0.0 

BO=C.C 

Ll= I SAVE(MM) 

no 120 I=1,K 

RM = RM+B ( I ) >"RY( I ) 

l=isavl ( n 

B( I )^B( I (STOILD/STDIL) ) 

120 BC=BC+B I I ) >^XBAR (L ) 

BO=XBAR (LI )-B0 
SSAF = RM>f^D( LI I 
122 RM= SQFTI ABS(RM)) 

SSDR=D( L 1 ) -SSAR 
FN=N-K-1 
SY=SSDR/EN 
DO 130 J=1,K 
L1=K^(J-1 ) +J 
L=ISAVE( J) 

125 Sb(J)= SOFT! A B SI ( R XI L 1 ) /D( L ) ) SY ) ) 

130 T( J) =BI J)/SB( J ) 

135 SY= SORT! ABSISY)) 

FK=K • 

SSAPM=SSAF/FK 
SSDRM=SSDF /FN 
F=SSAPM/SSDRM 
ANS(1)=BC 
ANSI 2) =PM 
ANSI ? )=SY 
ANSI^) =SSAP 
Ah SI 5) =FK 
ANSI 6 ) = SSAPM 
ANSI 7) = ssrp. 

ANSI 8) =FN 
ANSI 9) =SSDFM 
ANSI lCI^=f 
RETURN 
END 
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